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The spectral fluctuations of quantum (or wave) systems with a chaotic classical (or ray) limit are 
mostly universal and faithful to random-matrix theory. Taking up ideas of Pechukas and Yukawa 
we show that equilibrium statistical mechanics for the fictitious gas of particles associated with 
the parametric motion of levels yields spectral fluctuations of the random-matrix type. Previously 
known clues to that goal are an appropriate equilibrium ensemble and a certain ergodicity of level 
dynamics. We here complete the reasoning by establishing a power law for the % dependence of the 
mean parametric separation of avoided level crossings. Due to that law universal spectral fluctuations 
emerge as average behavior of a family of quantum dynamics drawn from a control parameter 
interval which becomes vanishingly small in the classical limit; the family thus corresponds to a 
single classical system. We also argue that classically integrable dynamics cannot produce universal 
' spectral fluctuations since their level dynamics resembles a nearly ideal Pechukas- Yukawa gas. 
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Most quantum systems whose classical limit is chaotic display universal spectral fluctuations [ij-Q. Such univer- 
sality was first observed in nuclear spectra, but is by now well supported by atomic and molecular spectroscopy. 
Somewhat surprisingly at first came support from microwave resonators and later even from elastic vibrations of 
crystals, i. e. classical waves for which the ray limit is chaotic. Such experimental as well as numerical findings suggest 
that discrete spectra of non- integrable quantum or classical waves display universal fluctuations. 
£-C) ■ The universal fluctuations in question can be roughly characterized as level repulsion. They are also met with in 
| spectra of random Hcrmitian or unitary matrices drawn from the Gaussian or circular ensembles of random-matrix 
theory Q pioneered by Wigner and Dyson. They contrast to the spectral fluctuations of integrable systems with 
two or more freedoms where we find level clustering, just as if the levels followed one another like the events in a 
Poissonian random process. 

The general rule just sketched is not without exceptions, some of which will concern us below. It is widely accepted, 
though, that hyperbolic classical dynamics can fail to being faithful to random-matrix theory (in their level-spacing 
> distribution as well as low-order correlation functions of the level density) only due to either quantum or Anderson 
localization or to slight perversity, like symmetries relevant for quantum spectra but not coming with classical constants 
of the motion Faithfulness to random-matrix theory even prevails for dynamics with mixed phase spaces, as 

long as islands of regular motion remain smaller than Planck cells. 

Theoretical attempts to explain the applicability of random-matrix theory to quantum and wave chaos have been 
based on treating the parametric dependence of spectra with the methods of statistical mechanics arL d, more 

recently, on extending ideas from the theory of disordered systems to dynamical systems • Building on the first 

d ' °f these strategies we here propose to put forth arguments which could lead to clarifying the issue. 

Parametric level dynamics was related by Pechukas and Yukawa |8) to the classical Hamiltonian motion of 
a one dimensional gas, the Pechukas- Yukawa gas (PYG) - Different PYG's arise Q for Hamiltonians depending on 
a parameter t as H(t) — Hq + tH\ or H(t) — HqCost + Hi sinr or H(t,r) = Hq + tH\ 2 — *)• The 

latter involves periodic kicks; the time evolution over one period is described by unitary Floquet operators of the 
form F = e _lff °e _1Tifl whose quasienergies ip m are defined by F\m,r) = e~ liPm \m, r). Note that we have taken the 
liberty to absorb a unit of time as well as Planck's constant in the Hamiltonians so as to make them as well as the 
perturbation strength r dimensionless. 

We find it convenient to base the following arguments on the PYG for periodically kicked systems (but we shall 
occasionally draw in illustrations for autonomous systems as well). The quasienergies (/3 m (r) with m = 1 . . . N appear 
as the coordinates of particles changing as the "time" r elapses. The role of conjugate momenta is taken up by 
Pm(T~) = (m,T\Hi\m,T). The Hamiltonian function comprises a kinetic term and a pair interaction [Q, 

n, 1 \ " 2 , 1 \ " \lmn\ 2 , 1 , 

m=l m^tn 
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As the only unusual feature of that Hamiltonian we encounter coupling strengths l mn which are themselves dynamical 
variables, related to the off-diagonal elements of the perturbation as l mn = i{m, r | J H r i|n,r)(e i( ^^^" ) - 1). If time 
reversal invariance holds for the quantum dynamics (for the other universality classes similar considerations apply), 
the initial Z mn (0) can be chosen to form a real antisymmetric matrix with N(N — l)/2 independent elements, l mn — 
—Inm = Imm an d the dynamics of the gas preserves that property. To extract Hamiltonian equations of motion for the 
fm,Pm, lmn wc need to employ the familiar Poisson brackets for canonical pairs {p m , <fn} = 6 mn ; the Z mra , on the other 
hand, Poisson commute with the ip m ,Pm and between themselves obey {l m n,kj} = ^(Smjlni+Snilmj—Snjlmi — fimilnj)- 
The latter Poisson brackets happen to be the ones for generators of rotations in iV-dimensional real vector spaces. 
The number of variables is D = 2TV + TV(TV — l)/2. Incidentally, if we were to solve the resulting Hamilton equations 
through a power series in the "time" t we would recover perturbation theory. 

Prior to appearing as a concise formulation of parametric level dynamics or perturbation theory, the PYG Hamil- 
tonian had been known to mathematicians and noted for being integrable, in spite of the non-trivial two-body 
interaction. It is easy to see that the PYG moves on an TV-torus in its Z?-dimensional phase space: In the eigen- 
representation of the perturbation, Hi\r) — u> r \r), the Floquet matrix takes the form F rs = (e~ lH °) rs e~ 1TU>! ' . The 
r-dependence thus lies in the TV exponentials e -I ™ s which define an TV-torus. Going from the _ffi-representation to 
the ^-representation amounts to a reparametrization of phase space by a canonical transformation which can deform 
but not topologically destroy the torus. If the TV eigenvalues cu r are incommensurate the TV-torus will be covered 
ergodically as the "time" r elapses. 

When applying equilibrium statistical mechanics to the PYG we must resist the temptation to employ the canonical 
ensemble p(<p,p, I) ex e - ' 371 ' since that ensemble allows the PYG to visit everywhere on the (D — l)-dimcnsional energy 
shell with uniform probability while only the much smaller TV-torus mentioned is accessible. It was nevertheless 
an important observation S that e~P u gives precisely the eigenphase density of random-matrix theory, p((p) = 
J d N pd N ( N ~ 1 ^ 2 l e~f 3u oc Ilm<n \e~ llfm — e _lv ™|, as is easily checked by doing the Gaussian integrals over the p's and 
Vs. The appropriate ensemble to use for the PYG is the generalized canonical ensemble 

p(<p,p,i) oc exp | - 2^ =1 /V^uj; ( 2 ) 

the accessible TV-torus is thus nailed down by fixing D — N suitable constants of the motion C u in the ensemble mean 
with the help of Lagrange parameters (3 U ; one of these C u should be the Hamiltonian H.. That ensemble has been 
shown ju| to have the distribution of level spacings as well as low-order correlation functions of the level density 
in common with random-matrix theory, to within corrections of order 1/TV. A good step towards explaining the 
universality of spectral fluctuations was thus taken, but several more steps are still ahead of us. 

As the phase-space trajectory of the fictitious gas winds around the TV-torus, the original Floquet matrix F{t) 
changes within a one-parameter family. It is that family rather than a single dynamical system (which has a fixed 
value of r) which is assigned random-matrix type spectral fluctuations. Indeed, ergodicity on the TV-torus implies 
that r averages of spectral characteristics like the spacing distribution P(S,t) equal ensemble averages, but the 
instantaneous form of P(S, r) is thus not specified. To improve the status of equilibrium statistical mechanics we 
inquire about the minimal time window At needed for time and ensemble averages to become equal. One might 
demand At — » oo. Such generous a provision would allow to accommodate and render weightless equilibration 
processes originating from initial conditons far out of equilibrium. As illustrated in Fig. |], equilibration does in fact 
take place in the spectrum of F(t) = e~ lH °e~ lrHl if Hq is integrable and chaos develops only as the perturbation is 
switched on; beyond a certain "relaxation time" the phase space regions of regular motion have shrunk to relatively 
negligible weight. If, on the other hand, Hq and Hi are both non-integrable (and from the same universality class) the 
initial state of the fictitious gas is already close to equilibrium, and then the window At in question need not be much 
larger than the collision time t co ii of the gas, i. e. the mean distance of avoided crossings for a pair of neighboring 
levels. 

We now propose to argue that T co n scales with the number of levels TV and thus, due to Weyl's law, with Planck's 
constant like a power, 

7-coii cx N- u « h fu , v > 0, (3) 

where / is the number of degrees of freedom. A T-average for the PYG over a window At oc t co u thus involves a 
family of Floquet operators F(t) which all y ield identical classical dynamics in the limit Ti oc TV~7 — > 0. That insight 
grew out of an idea put forward in Ref . [|ll| . 

The following estimates support the power law (^|) and reveal the exponent v as non-universal. First, consider a 
Hamiltonian H(t) = Hq cost + H\ sinT w Hq + tH\ and let Hq and H\ be independent random TV x TV matrices 
drawn from the Gaussian orthogonal ensemble. Both should have zero mean, H oi j = Hnj = 0, and the same variance, 
\Hoij\ 2 = \Huj\ 2 = 1/TV. The mean level spacing then is A = TS^+i — Ei = O(jj). The level velocity vanishes in the 
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ensemble mean, Hm — 0, and has the variance p 2 = H 2 U = O(jj). The estimate r co n ~ A/p ccJV 2 thus yields the 
exponent v — \. The same value results for Floquet operators F(t) = e _lH °e~ 1THl with Ho, Hi random as before. 

Another example is provided by the kicked top Q of Fig. [j] for which the operators Hq and Hi are functions of an 
angular momentum J such that J 2 = is conserved; for fixed j the dimension of the Floquet matrix is N = 2j+l. 

If we choose the perturbation as a torsion, Hi = J 2 /(2j + 1), the mean level velocity is V = -k Yli=i Hw = ^TriJi = 
Y^(l + O(jj)); this gives a common drift of all levels irrelevant for collisions. A typical relative velocity p may be 
found from the variance P 2 = 5"^— i (Hm — v) 2 , where the matrix elements are meant in the eigenrepresentation of 
-F(t); but in that representation the perturbation Hi = Hi —v will, for chaotic dynamics, look like a full matrix with 
E l Hm » and £. » £ E« = ^ Tr ^? = + °(^)); we conclude p cx iVi The mean level spacing ^ 
yields the collision time cx N~i and the exponent v = |; we have confirmed the latter value by numerically following 
level dynamics for 10 < j < 160. 

Somewhat daring but interesting are the following estimates for some classes of Hamiltonians with / freedoms. It 
is not untypical to find H represented by a sparse banded N x N matrix with a zebra-like structure of non- vanishing 
elements. An example for such a Hamiltonian is given by a finite cubic lattice with TV = M' sites and nearest 
and next-to-nearest neighbor interactions. Non-vanishing matrix elements then only exist in approximate distances 
M°,M 1 , . . . ,M^ _1 from the diagonal. The overall bandwidth is thus AP^ 1 = JVu -1 )//. For similar Hamiltonians 
we have numerically studied the collision time, choosing H(t) = H n cost + Hi sinr as functions of SU(3) generators, 
like the Lipkin-Meshkov-Glick Hamiltonian known from the nuclear shell model Jig ]. Such SU(3) dynamics can also 
arise for collections of identical three-level atoms collectively interacting with (nearly) resonant modes of a microwave 
or optical cavity |l6| ], and such realizations could even access the (semi) classical limit as the number of atoms grows 
large. That limit can, for a given Hamiltonian, lead to classical motion with either / = 2 or / = 3, depending on 
the representation the initial state belongs to. We have chosen a Hamiltonian yielding global chaos in both cases and 
numerically followed level dynamics over two decades of N (see Fig. ^). The collision time exponents v were measured 
using T co u = A/p and were found as 0.737 in the two-freedom case and 0.677 for / = 3. These numerical results lead 
us to conjecture that the collision time for this class of Hamiltonians has the exponent Vf — The limit / 3> 1 

corresponds to the above estimate v = i for full random matrices. 

We now turn to integrable dynamics and their PYG's. The Hamiltonian ([l]) still applies. Only the initial conditions 
for <fr,p, I, i.e. the matrices Hq, Hi, tell the PYG that the classical variant of F(t) is integrable. We propose to show 
that the PYG then tends to behave like a near ideal gas. Assuming / freedoms and the Hamiltonian expressed in 
terms of action variables Ji, . . . ,It the semiclassical Einstein-Brillouin-Keller (EBK) approximation for the spectrum 
results from letting each action run through integer or half integer multiples of Planck's constant, li = H(rii + Mi/4), 
with rii = 0, 1, . . . and fii the Maslov index. A multiplet of levels arises as one of the / quantum numbers runs while 
all others are fixed. EBK levels of different multiplets have no reason not to cross as a control parameter r is varied. 
The exact levels will also cross, provided there are / — 1 independent operators commuting mutually and with the 
Hamiltonian. Assuming such quantum integrability, we may treat intra- multiplet level dynamics for / = 1. 

To elucidate the level dynamics for a single-freedom system we may first look at a non-symmetric double-well 
potential with a finite barrier separating the wells; there will be a single eight-shaped separatrix in the phase space, 
with the energy of the top of the barrier. For energies above the (r-dependent) barrier, all classical orbits lie outside 
the separatrix and the function H{I , r) is single valued and monotonic. Variation of r cannot bring about crossings 
of quantum levels, exact or semiclassical, which stay above the barrier. On the other hand, for energies below the 
barrier, separate classical orbits exist within each loop of the separatrix and for the quantum energy cigenfunctions 
are similarly localized. The function H{I) now has two branches each giving a different action for a given energy. 
Neighboring EBK levels whose eigenfunctions live in different wells can be steered through a crossing by varying 
r; only tunneling corrections turn such crossings into avoided ones. The closest-approach spacing in an avoided 
crossing is determined by the (modulus of the imaginary) action AI — j dxy // 2m(E — V{x)) across the barrier as 
~ e _ ' A/ '/ ?i ; that spacing may be unresolvably small for \AI\/fi 3> 1. In the language of the fictitious gas we could 
speak of extremely weak repulsion or near ideal-gas behavior. A few EBK crossings can become, by tunneling, strongly 
avoided, with closest approaches of the order of the mean spacing; these appear for levels near the top of the barrier 
where AI/h is of order unity. 

Fig. displays the level dynamics of a single-freedom system with finite N and illustrates the foregoing discussion. 
For t below a certain critical value nothing like a crossing occurs since H{I,t) remains monotonic and single valued 
in /. For r > r c very narrowly avoided crossings appear which on the scale of the plot look like crossings. In the 
neighborhood of a curve E = £ , sep (r) giving the energy of the separatrix, not drawn in the graph but clearly visible 
nonetheless, we see strongly avoided crossings. Inasmuch as a finite fraction of all N levels can hit the curve E sep (T) we 
may argue that the number of strongly avoided crossings is linear in N while the number of unresolved anticrossings 
is ~ A^ 2 . Of course, the behavior just sketched pertains to Hamiltonians with a single separatrix. We could easily 
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construct a Hamiltonian with so many separatrices and consequently so many strongly avoided crossings that its level 
dynamics for fixed matrix dimension N would look rather like what we usually find for chaotic systems; but such a 
Hamiltonian would at once be revealed as an impostor by decreasing the effective value of h, i. e. going to matrix 
representations of larger dimension N; the very mechanism making more or less all crossings well avoided for some 
sufficiently small N would cast them into negligible minority as N — > oo or, equivalently, % — > 0. 

We wish to sum up our findings. Under conditions of chaos the mean parametric distance of avoided level crossings 
goes to zero in the semiclassical limit like a power of h. Consequently, an average over a classically vanishing control 
parameter interval reveals universal spectral fluctuations. For classically integrable dynamics, on the other hand, level 
repulsion is non-generic. 

Open to further study is the question as to which, if any, equilibrium ensembles of the PYG apply to classically 
integrable dynamics and to chaotic dynamics not faithful to random-matrix theory, like those with (Anderson) local- 
ization, quantum symmetries without classical counterpart HQ or otherwise intermediate statistics Jl8| . To specify 
such ensembles we would need a better understanding of the (infinity of) constants of the motion of the PYG. 

We have enjoyed discussions with Joachim Weber, Yan Fyodorov, Jon Keating, Hans-Jurgen Stockmann O], and 
Martin Zirnbauer, as well as support by the Sonderforschungsbereich "Unordnung und Grofie Fluktutionen" der 
Deutschen Forschungsgemeinschaft and Polish KBN Grant No 2 P03B 023 17. 
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FIG. 1. Quasienergy levels tp of the kicked top with the Floquet operator F(t) = e ~ l,3Jy e~ laJx e~'" I 3+ Tj * 2 for j = 15, 
a = 0.907, (3 — 1.072. As r grows from to ~ 5 the spectrum "relaxes" from level clustering to level repulsion, while the 
dynamics goes from integrable to chaotic. 
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FIG. 2. A^-dependence of the level collision "time" k Af " for a Hamiltonian H(t) = Hq cost + i/isinr with Ho, Hi 
polynomials in SU(3) generators. Such an H(t) can correspond to either 2 or 3 classical freedoms. 
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FIG. 3. Level dynamics of the integrable spin Hamiltonian H(t) — (J z + 0AJ x ) cost + (jf^jJz ~ 0-5J z — ^f-) sinr with 
j — 15. Strongly avoided crossings are rare exceptions. Right box is a blow-up of left. 
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